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ABSTRACT 

We have performed a series of chemodynamical simulations of barred disc galaxies. 
Our goal is to determine the physical processes responsible for the increase in the 
central gas-phase metallicity and of the central star formation rate (SFR) observed 
(~| ' in the Sloan Digital Sky Survey (SDSS). All simulations start with an axisymmetric 

Qlh! distribution of stars and gas, embedded into a spherical dark matter halo. We define 

a 2kpc diameter central aperture to approximate the integrated spectroscopic fibre 
measurements from the SDSS. The chemical evolution observed within this central 
, region depends critically upon the relative size of the bar and the aperture, which 

^ ' evolves strongly with time. At t 0.5 Gyr, a strong bar forms via a disk instability, 

whose length is considerably longer than the 2kpc aperture. The stars and gas lose 
angular momentum and follow elongated orbits that cause an intense mixing of the 
gas between the central region and its surroundings. During the next 1.5 Gyr, the 
orbits of the stars inside the bar do not evolve much, but the orbits of the gas contract 
significantly until the entire gas bar is contained in the 2kpc aperture, resulting in 
' a net flux of gas into the central region. During this period, the metallicity in the 

\^ , central region increases steadily, and this enrichment is dominated by metal-rich gas 

that is flowing into the central region. The main result of this work is therefore that 
the observed enrichment in the centres of barred galaxies is not dominated by in-situ 
. enrichment by stars formed in the centre. Rather, star formation occurs along the full 

length of the bar, much of which occurs initially outside the 2 kpc aperture. About 50% 
of the metals that end up in the central region originate from this extended bar-long 
star formation, but flow into the central region due to loss of angular momentum. The 
K> ' effect is less significant for iron because the delay for the onset of Type-la SNe leaves 

;_j , less time for mixing. Still, there is a significant increase in [Fe/H] before the central 

' stars contribute to in-situ enrichment. Eventually, as the orbits of the gas inside the 

bar contract, they fall completely inside the central region, and only then the central 
region can be regarded as a closed-box system. However, by that time, most of the 
metal-enrichment in the central region has already taken place. We conclude that 
there is no direct connection between central SFR and central metallicity. The central 
metallicity does not originate exclusively from central stars. Instead, the global SFR 
(especially along the bar) and the large-scale flow of enriched gas play a major role. 

Key words: galaxies: evolution — galaxies: spirals — ISM: abundances — stars: 
formation 
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1 INTRODUCTION 

1.1 Evolution of Barred Galaxies 

Bars are a very common feature in spiral galaxies. Ob- 
servational determinations of the bar fraction in nearby 
galaxies, and up to redshift z ~ 1, systematically 



produce values in t h e range 20% to mor e than 60% 
jEskridge et al.1 l2000l: lElmegreen et al.1 12004 iJogee at al.l 
2004'; 'Nair & Abraham' '2010"; 'Martmcz & Muriel' '201 1|; 
Giordano et al. 2011; Masters ct al. 2011; Lee et al. 2012ai). 
The formation of these bars can be trigg ered by merger s 
or tidal interactions with other galaxies (|Noguchil Il987l ') , 
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or can happen in isolated gal axies as a result of dynami- 
cal instability (^Sellwoodlll981^ ■ Analytical studies and nu- 
merical simulations have shown that the presence of a bar 
can greatly affect t he evolution of the host galaxy. Ana - 
lytical calculations (|Lvnden-Belllll979l : lAthanassoulall2003l ) 
showed that the bar produces a gravitational torque that 
transports angular momentum outward and matter inward. 
This was confirm e d by t he ea rly simulations of Roberts et al.l 



1981 ) 



19911 ). 



l| 19791 ) . ISellwoodl (|l98ll ). and Ivan Albada fc Roberts 

In the more rec e nt simulatio ns of|pfcnnigcr fc Friodli ( 

iFriedli fc Be"n3 (|l993l . Il995') and lFriedli et all (|lS?94 ). a bar 
forms by dynamical instability. That bar channels gas into 
the center of the galaxy, resulting in a flattening of the 
initial metallicity gradient. Eventually, the gas flow causes 
a central starburst, and the resulting metal enrichment 
steepens the metallicity gradients in the central regions. 
Several autho rs studied the dynamic s of the gas flowing 



several autho rs studied the dynamic s ot the gas nowmg 
along the bar JCombes fc"G erin 198 51: lAthanassoulalll992b|: 
Elmegreen' 'l993; Fu xi llS)99l: iMacieiews ki et all 



Combes 

2002; Re gan fc Toubcn 2004 ; Baba et al.1 [20101. Their sim- 
ulations showed that the gas retains some angular mo- 
mentum, preventing it from falling directly into the cen- 
ter of the galaxy. Instead, the gas tends to form elongated 
orbits inside the stellar bar. This result is in agreement 
with analytical studie s of the properties of orbits in ax- 
isymmetric potentials (|Athanassoulain"992al ) . If gas does fall 
into the center of the galaxy, it might acc rete onto a cen- 
tral black hole and fuel a cent r al AGN |Shlosman et al.l 
1989'; 'Shlosman fc Noguchi' fl993 : iHeller fc Shlosmanlll993 : 
Combes. . 2003.; .Jogeg ,2006). Or, it might lead to the f orma- 



tion of a central bulge ( Kormendv fc Kennicutt]|2004 ) 

All these analytical studies and numerical simulations 
provide a basic scenario for the evolution of barred galax- 
ies. Once the stellar bar forms, gas starts flowing along the 
bar, into the central regions. This inflow of metal-poor gas 
tends to flatten the metallicity gradient. The star formation 
rate in the central regions increases, either steadily or in the 
form of a starburst. This will lead eventually to an increase 
in central metallicity. Finally, the gas that is not converted 
to stars might end up fueling a central AGN. We would 
therefore expect barred galaxies to have higher central SFR, 
higher central metallicities, higher bulge fraction, and more 
intense AGN activity than unbarred galaxies, and possibly a 
flatter metallicity gradient, depending at what stage of their 
evolution they are being observed. This basic scenario is in- 
deed supported by many observations. Some barred galax- 
ies have shallower metallicity gra dients than unbarred ones 
l|Pagel et al.ll 19791 : iMartin fc r"wi,1994 1 , suggesting gas flow 
along the bar. Several observations reveal enhance d star for- 
mation in the central regions of barred galaxies ('Heckman' 
1980; Hawardcn ct al . 1986 ; Dcvcrcux 1987; Arscnault 1989; 
Huang et al.lll996l: IHo et ar.lll997l: iMartinet fc Friedlilll99^: 



Hunt fc Malkanlll999l: lEmsellem et al.ll200ll: iKnapen et al 



2002 



200S 



Laurikainen et al jl2004l : iJogee et al.ll2005ir iHunt et al. 
Bang fc Annll2009l ). Also, [Masters et all l|2012h found 



an anti-correlation between the presence of a bar and atomic 
gas content, which could be interpreted as a depletion of 
gas in barred galaxies, possibly res ulting from enhanced 



star formation. [Masters et al.l |20l3) found that over half 
of red, bulge-dominated disc galaxies possess a bar, while 
blue galaxies show no evidence for a bar or a bulge. The is- 
sue of AGN fueling is less clear. Several observational stud- 



ies find either a larger bar fraction in AGN-host galaxies 
than non-AGN ones ( fArsenault 1 989 : iKnapen et aU I2OO0I : 
iLaine et al.ll2002l ) or a larger AGN fraction in b arred galax- 
ies than unbarred ones (jCoelho fc Gadottill2011 ), other stud- 
ies do not find any significant difference jMoles et al.l 19951: 
McLoed fc Riekelll9 95': 'Mulchaev fc Regan" 'l99 7l: IHo et al 



1997 : Xaurikainen et al.,,2004 : ,Hao et al..,2009. : .Bang fc Ann 
2009l : lLee et al.ll2012bl '). 



1.2 Recent Observations 

While many observations support the basic scenario of bar 
evolution described above, several recent studies of barred 
galaxies using large samples reveal a more complex pic- 
ture. A study of 294 galaxies with s trong bars fr o m the 
SJoan Digital Sky Survey (SDSS) by lEUison et all (|201ll ) 
showed that barred galaxies of stellar masses M. > 10" M0 
have enhanced central SFR and metallicities, which puts 
the m as outliers on the fun damental m etallicity relation 
(.Ellison et al.ll2008l : iMannucci et al.. 201ol : iLara-Lopez et all 
l2010al lb^. Galaxies with stellar masses M« < 10^" also 
show an increase in central metallicity, but without a corre- 
sponding increase in central SFR. One possible explanation 
is that star formation in the center of low-mass barred galax- 
ies, which is presumably responsible for the observed higher 
metallicities, has now ceased, while stars are still forming in 

the center of high-mass galaxies. 

More recently, Wang et^l] l|2012l ') performed a detailed 
study of a sample of 3757 galaxies from SDSS, including 
1555 barred galaxies. They found that the presence of a 
bar does not automatically imply higher SFR. Galaxies with 
weak bars (ellipticities Cbar < 0.5) do not show enhanced 
central star formation, only galaxies with strong bars do. 
These strong bars are found predominantl y in galaxies with 
high stellar masses M, > l-3x 10^°Mm iNair fc AbrahamI 
l20ld : lEUison et aT]|201ll : IWang et al.ll2012l '). This implies a 
relation between mass and S FR, which could explain the 
results of lEllison et al.l (|201lh . B ut interestingly, wh ile only 
galaxies with strong bars in the IWang et al.l 1 20121 ) sample 
have enhanced central SFR, not all of them do. Instead, 
some of these galaxies have a central SFR that is lower than 
the mean by a factor of 10. This suggests that bars (or some 
other process associated with their presence) might, under 
some circumstances, quench star formation in the central 
regi ons. 

lOh et aU (|2012l ) studied the effect of bars on central star 
formation and AGN activity, using a sample of 6658 late- 
type galaxies from SDSS, including 2442 galaxies with bars. 
They find that the presence of a bar significantly increases 
the central SFR in red galaxies, but not in blue ones, where 
the central SFR is large whether or not a bar is present. This 
could result from the fact that redder galaxies tend to have 
longer bars, which according to simulations leads to larger 
effects. Or, as the authors suggest, central star formation is 
naturally reduced in redder galaxies because less gas is avail- 
able, making the effect of infalling gas more important than 
in bluer galaxies. While bar effects on the central SFR are 
important for red galaxies but not blue ones, the opposite 
is true for AGN activity. Bar effects on AGNs are mostly 
found in blue galaxies, even though most AGNs are hosted 
by red galaxies. The authors suggest that the presence of 
a large central mass concentration can weaken the bar, re- 
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ducing its effects. Another possibility is that red galaxies 
hosting AGNs are in the post-starburst phase. Then most 
of the gas has already been co nsumed by star fo rmation, 
leaving no gas to feed the AGN (|Sheth et al.ll2005l 'l. 

These recent studies reveal that the evolution of barred 
galaxies is probably more complex than the basic scenario 
suggests. Galaxies of different masses, or with bars of dif- 
ferent strength, can have drastically different star forma- 
tion histories. As a result, there is no simple correlation 
between the presence of bars, enhanced central SFR, and 
enhanced central metallicity. To fully understand the evolu- 
tion of these galaxies, we need to understand the interplay 
between the gravitational dynamics responsible for the for- 
mation of the bar and the large-scale gas flows in the galaxy, 
the hydrodynamics responsible for the formation of stars and 
the chemical enrichment of the ISM, the feedback processes 
resulting from stellar evolution and AGN activity, and their 
possible role in dispersing chemical elements in the galaxy, 
regulating star formation, and possibly quenching it. 



1.3 Objectives 

To address this complex problem, we first focus, in this pa- 
per, on the connection between enhance d central star forma - 
tion and enhanced central metallicity f Elliso n et al.ll201ll ). 
This is a first step toward understanding the evolution of 
barred galaxies and building a scenario which is consistent 
with all recent observations. Previous hydrodynamical sim- 
ulations of barred galaxies showed that gas fiow along the 
bar usually results in increased st ar formation and increased 
metallicity in the central regions (iPfenniger fc Friedlil Il99ll : 



iFriedli fc Benj|l993l . [l995l : lFriedfi et al.lll994l ). butlhat does 
not imply that one is the cause, or partial cause, of the other. 
Metals that end up in the central regions might have been 
produced elsewhere in the galaxy, while metals produced 
in the central regions might eventually leave that region. 
These specific questions were not addressed in previous sim- 
ulations. 

We have performed a series of chemodynamical simu- 
lations of isolated barred galaxies. We focus on the evolu- 
tion of the metal abundance in the central region. We define 
this central region as the central 2kpc, because it roughly 
matches the size of the fibre region in the SDSS. From now 
on, we will use the term "central," to refers to the 2kpc 
aperture, except in section 3.6 where we consider different 
aperture sizes. We identify all physical processes responsi- 
ble for modifying the abundance of individual chemical el- 
ements in the central region, and we determine the relative 
importance of these various processes. Our objective is to 
determine if enhanced central star formation is the primary 
cause for the enhanced central metallicity, or if the latter has 
a different origin. This is an essential first step in assessing 
the validity or shortcomings of the basic scenario. 

The remainder of this paper is organized as follows: In 
§2, we describe the numerical algorithm used for this study. 
Results are presented in §3, where we first discuss the global 
evolution of the galaxy, and then focus on the central region. 
The implications of the results are discussed in §4. Summary 
and conclusions are presented in §5. 



2 THE SIMULATIONS 

2.1 The Numerical Algorithm 

We us ed the chemodynamic a l numerical algorithm 
GCD+ (iKawa-ta fc GibsonI l2003l : iRahimi fc Kawatal I2OI2I : 
iKawata et all |2013| ). GCD-I- is a three-dimensional tree 
N-body/SPH code that incorporates self-gravity, hydro- 
dynamics, radiative cooling, star formation, supernova 
feedback, metal enrichment, and metal diffusion. Gas is 
converted into stars in regions where the number density 
exceeds a certain threshold nth and where the velocity 
field is convergent , follow ing the Schmidt law described in 
iKawata fc GibsonI l|2003l ): 

dp, _ _dpg_ _ CtPg 

dt dt t„ 



(1) 



where pg and p» are the gas and stellar mass density, re- 
spectively, C» is a dimensionless SFR efficiency, and tg is 
the dynamical time. 

We assume that the the stel la r ma sses are dis- 
tributed according to the ISalpeteJ l| 19551 ) initial mass 
function (IMF). Ch e mical enrichment by both Type II 
i Wooslev fc Weaverl [l995l ) and Type la supernovae 
l lwamoto et al .' 1999') and mass loss from intermediate- 
mass stars l|van don Hock fc Gro cncwcgcn 1997) are taken 
into account . We adopt th e Typ e la supernovae rate 
suggested bv lKobavashi et al.l (|2000l ). We assume that each 
supernova produces an amount of thermal energy -Esn, 
and that stellar winds from massive stars (M > 30A'/q) 
also produce thermal energy, at a rate -Esw. The canonical 
energy of a supernova is 10^^ erg. However, according 
to the high-resoluti on one-dimensional simulations of 
iThornton et all l| 19981 ). 90% of the initial supernova energy 
is lost in radiation in its early expansion phase, and does not 
contribute to feedback. Accordingly, we set Esn = lO^^erg 
per supernova for our fiducial model. Notice that an initial 
supernova energy has not been established quantitatively 
yet. For this reason, we also explore the effect of a higher 
supernova energy, _Esn = 10^^ erg, in one of our simulations. 
We set -Esw = lO'^^ergs"^ for our fiducial model, but 
also consider the effect of a stronger stellar wind with 
^5sw = l O^^ergs"'^, to cover t he range sugge sted by s everal 
authors (IWeaver et all Il977l : IShuU 1980,; iGibsonI Il994l : 
lOev fc Massevlll994) . The metal diffusion is modeled with 
the method described in iGreif et al. I l|2009l ). We use a star 
formation density threshold, nth 
formation efficiency, C, = 0.02. 



0.1cm and a star 



2.2 Initial Conditions 

The technique use d for generat i ng ini tial conditions is de- 
scribed in detail in lGrand et al.l (|2012l ). We use particles to 
represent the stellar and gaseous components of the galaxy. 
The stella r and gas discs are s et up using the method de- 
scribed in lSpringel et al] (|2005l ). The stellar disc follows an 
exponential surface density profile: 



Pd,sta 



Attz /?2 

^" "^d, stars -"'d, stars 



-R/ -Rd.st, 



' sech 



(2) 



where R, z, Md, -Rd and the radial coordinate (dis- 

tance from the z-axis), vertical coordinate, disc mass, scale 
length, and scale height, respectively. The gaseous disc has a 
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Table 1. Numerical Parameters of the Initial Conditions 



Table 2. Numerical Parameters of the Simulations 



Disc 


Md [Mq] i?d [kpc] 


Zd [kpc] 


stars 


4 X 10^0 2.50 


0.35 


gas 


1 X 10^0 4.00 





same radial dependence as the stellar disc, but with a differ- 
ent scale length 7?d,gas, and its initial vertical distribution 
is not given by a sech z profile. Instead, it is adjusted by 
imposing that the gas is initially in hydrostatic equilibrium. 
The values of the parameters are given in Table [T] These pa- 
rameters were c hosen in order to produce a galaxy similar to 
the Milky Way (|011ing fc Merrifieldll200ll : iRobin et al.ll2003l : 
lMcMilIan|[20Tll 'l. In this first study, we did not include gas- 
infall from the surrounding intergalactic medium, in order 
to focus on the effect of gas accretion in the disc. 

We impose initial radial abundance profiles to both stel- 
lar and gaseous components. The iron abundance profiles for 
the stars and gas are given by 



[Fe/H](i?) = 0.2-0.05 



R 
1 kpc 



The a-element abundances are given by 
-0.16[Fe/H], stars; 



[a/Fe](i?) 



.0, 



gas. 



(3) 



(4) 



At each radius, we add to the values given by equations ((3]) 
and Q a gaussian scatter with a dispersion of 0.02 dex. For 
the star particles, we use the value of [Fe/H] to randomly 
assing them an age, following the age-metallicity relation 
[Fe/H] = —0.04 X age (Gyr) with a 0.02 dex scatter, and 
assuming a constant SFR in the — 10 Gyr age range. 

We do not use particles to represent the dark mat- 
ter halo. Instead, we impose a fixed, non-evolving spherical 
gravitational potential, in order to focus on the baryon dy- 
namics and save computational time The dark matter halo 
density follows a truncated NFW profile. 



(r/r,)(l+r/r.)2^ 

l|Navarro et al.l Il997l : iRodionov et all |2009| ) where r is the 
distance from the center. The scale length is given by 
r20o/c, where c is the concentration parameter, and r2oo is 
the virial radius, inside which the mean density of the dark 
matter is equal to 200 times the the critical pcrit. The charac- 
teristic density pc is adjusted such that the mass M200 of the 
dark matter inside the virial radius is equal to 2 x IO'^^Mq. 
The truncation factor exp(— r^/r2oo) is introduced in our ini- 
tial condition generator for a live halo simulation. Although 
we use a static dark matter halo in this paper, we still use 
the profile given by equation ((5)l . Note that adding this trun- 
cation factor to the standard NFW profile leads to very little 
change in the central region, which we focus on in this paper. 



2.3 Runs and Parameters 

We performed a total of 4 simulations. The parameters of the 
simulations are listed in Table O We first performed a fidu- 
cial run, or Run F, with a medium resolution (120,000 star 
particles and 30,000 gas particles), concentration parameter 



Run 


-^stars 




c 


-E^SN [erg] 


-Bsw [ergs l] 


F 


120,000 


30,000 


10 


1050 


10^6 


HR 


360,000 


90,000 


10 


1050 


10^6 


NoBar 


120,000 


30,000 


20 


1050 


103*5 


Feed 


120,000 


30,000 


10 


1051 


1037 



nova, and stellar wind energy input -Bsw = 10 ergs"'^. We 
then performed a series of simulations which all differ from 
the fiducial run in one particular aspect: Run HR is a high- 
resolution run, with three times more particles. In Run No- 
Bar, we increased the concentration parameter of the dark 
matter halo, to prevent the formation of a bar. In Run Feed, 
we increased amount of stellar wind and SNe feedback by a 
factor of 10. All simulations end at t = 2 Gyr. 



3 RESULTS 

3.1 Global Properties 

3.1.1 The Formation and Evolution of the Bar and Spiral 
Pattern 

Figure [l] shows the distribution of stars and gas particles 
at six different times, for Run F. At t = 0.30 Gyr, we see a 
spiral pattern with multiple spiral arms, and a central bar is 
starting to form. At t = 0.50 Gyr, the central bar is fully de- 
veloped, and the spiral pattern is reduced to two main spiral 
arms. The stellar and gaseous components of the bar have 
a comparable length at that time. By t = 1 Gyr, the bar is 
still present, but the spiral pattern starts to disappear due to 
the dynamical heating by th e bar and spiral arms themselves 
(|Carlberg fc SellwoodiligSSl ). At t = 1.5 Gyr, the spiral pat- 
tern is hardly visible in the stellar component, while in the 
gas component we see several weak spiral arms. At the end 
of the simulation, the spiral pattern is gone but there is still 
a strong bar in the stellar component. Between t = 0.50 Gyr 
and the end of the simulation, the length and shape of the 
stellar bar remain essentially unchanged, while the gaseous 
component of the bar steadily contracts. At t = 2.00 Gyr, 
the stellar and gaseous components of the bar are 7 kpc and 
1 kpc long, respectively. 

We plot in Figure [2] the velocity field of the gas in the 
central region. On each panel, the red circle has a diameter 
of 2 kpc, and shows the central region. At t = 0.5 Gyr, the 
gas inside the bar moves along an elongate d, elliptical orbit. 
This i s consistent with the simulations of iRegan fc TeubenI 
(|2004h . This orbit intersects the edge of the 2 kpc aperture, 
hence the gas is flowing across the boundary, both moving 
in and moving out of the simulated flbre. This elongated 
vortex contracts with time (see Fig. [1} . At t = 1 Gyr, the 
gas bar is barely contained within the 2 kpc aperture, and at 
t — 1.5 Gyr, it is entirely contained inside the central region. 



3.1.2 Star Formation and Chemical Evolution of the 
Galaxy 

Figure |3] shows the time-evolution of global properties of the 
galaxy. The top left panel shows the SFR. The total galactic 



© XXX RAS, MNRAS 000,[THl8] 



Star Formution and Metallicity in Barred Spiral Galaxies 5 





-10-5 5 10-10-5 5 10 -10-5 5 10-10-5 5 10 
Ckpc) (kpc) (kpc) (kpc) 

Figure 1. Distribution of star particles (left panel in each pair) and gas particles (right panel in each pair) for Run F at six different 
times, as indicated. Each panel is 25 kpc X 25 kpc in size. 



SFR is initially constant at around AMq yr~^ . Then, at t = 
0.5 Gyr, the bar forms and starts chanelling gas in the central 
region, resulting in additional star formation. The SFR keeps 
increasing, reaching 6Mq yr~^ at t = 1 Gyr, and then starts 
decreasing. The decrease in SFR after t = 1 Gyr is not due 
to bar disruption, since the bar remains strong until the end 
of the simulation. It is simply due to gas exhaustion: star 
formation reduces the amount of gas available, and this gas 
is not replenished since we do not include gas infall from the 
intergalactic medium in our simulations. The bottom left 
panel shows the total hydrogen mass Mh in the gas phas^. 
It decreases almost linearly, with a slight change of slope at 
t = 1.1 Gyr. By the end of the simulations, it is down by 
50%, and correspondingly, the SFR is also down by about 
50% 

The right panels show the evolution of [O/H] and 
[Fe/H]. The evolution of [O/H] is closely related to the evo- 
lution of the SFR. The oxygen abundance rapidly increases 
during the period of high star formation, with [O /H] increas- 

^ Unless specified otherwise, elemental masses, abundances, and 
abundances ratios always refer to the elements present in the gas 
phase, not the stellar phase. Also, we do not distinguish between 
molecular, atomic, and ion species. Hence, Afjj refers to the total 
mass in the form of HI, HII, and H2. 



ing by 0.21 dex during the first f ~ 1.1 Gyr. Then, as the SFR 
starts to decrease, there is a significant change of slope, and 
[O/H] only increases by 0.04 dex during the last 0.9 Gyr of 
the simulation. The iron abundance increases more steadily, 
and at the end of the simulation [Fe/H] is increasing faster 
than [O/H]. This is caused by the different lifetimes of Type 
II and Type la SNe progenitors. Type II SNe have short pro- 
genitor lifetimes, between 5 and 10 million years. Since these 
SNe are the main producer of oxygen, we expect to see a 
tight correlation between the histories of star formation and 
oxygen enrichment. Iron is produced mostly by Type la SNe, 
which have long progenitor lifetimes, typically between 1 and 
13 Gyr. As a result, a significant fraction of iron is produced 
long after the SFR has reached it s peak. These resu l ts are 
consistent with the early results of lKawata fc GibsonI (|2003l ) 
(see their Fig. 2). 

Figure 2] shows the evolution of the radial profiles of 
[O/H] and [Fe/H], where R is the distance measured from 
the rotation axis of the galaxy. We calculated these profiles 
first by averaging over all angles, and then by only consid- 
ering the gas located in a cylinder of radius 1 kpc centered 
on the bar (solid and dotted lines, respectively). The black 
lines show the linear gradient present in the initial condi- 
tions. Once the bar forms, the fiow of gas along the bar tends 
to fiatten the abundance gradient. In particular, [O/H] and 
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Figure 2. Gas velocity in the central region for Run F at four different times, as indicated. The red circle shows the central region, and 
has a diameter of 2kpc. For clarity, only 1/4 of the particles are plotted. 



[Fe/H] in the central region initially go down, as metal-poor 
gas is channeled into the central region by th e baj. The same 
phenomenon was se e n in the si mulations of FriedU fc Bend 
(|l993l ). lFriedh et all (|l994l '). and lFriedli fc Bend (|l995l '). This 
is partly compensated by metals produced by Type II SNe. 
Since these SNe produce more oxygen than iron, we see a 
larger drop in the central [Fe/H]. The SFR increases with 
time, resulting in more SNe, and eventually the oxygen and 
iron abundances increase with time at all radii. The profiles 
tend to be flatter in the center when we consider only the 
gas located along the direction of the bar. The length of 
the plateaux seen at small radii in Figure U roughly corre- 
spond to the extent of the gaseous component of the bar, 
which is indicate by the vertical lines (see also Fig. [9] be- 
low. Note that at = 0.1 Gyr, the bar has not formed yet). 
Gas flow and strong mixing of gas inside the bar tends to 
erase any signiflcant abundance gradient. The plateaux get 
shorter with time simply because the gaseous component of 
the bar contracts. 



3.2 Inside the Central Region 

We deflne the central region as a cylinder of radius 1 kpc, 
parallel to the rotation axis of the galaxy, but centered on 
the center of mass of the stellar component. We chose this 
particular size because it corresponds to the typical flbre 
values sampled by the SDSS. 



3.2.1 Star Formation and Metal Abundances 

Figure [5] shows the evolution of the SFR, hydrogen, oxygen, 
and iron mass, and [O/H], [Fe/H], and [0/Fe] inside the 
central region. The dotted lines show the global evolution, 
copied from Figure [S] The chemical evolution of the central 
region differs signiflcantly from the global evolution of the 
whole galaxy. Initially, the amount of hydrogen in the central 
region is small, and correspondingly the SFR is also small. 
As the bar starts channeling gas toward the center, both Mh 
and the SFR increase, especially after t = 0.5 Gyr, when the 
bar is the strongest. At t ~ 1 Gyr, the gaseous component of 
the bar has contracted sufficiently to be entirely contained in 
the central region. The inflow of hydrogen stops, and both 
Mh and the SFR steadily decrease afterward. After t ~ 
1.2 Gyr, star formation takes place almost exclusively in the 
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Figure 3. Time-evolution of global quantities, for Run F. Top left panel: Star formation rate; bottom left panel; total hydrogen mass in 
the gas phase; right panels: [O/H] and [Fe/H] in the gas phase. 
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Figure 4. Radial profiles of [O/H] and [Fe/H] in the gas phase 
at various epochs, for Run F. Top panel: [Fe/H]; bottom panel: 
[O/H]. Colors represent different epochs, as indicated. Solid lines: 
azimuthally-averaged profiles; dotted lines: profiles calculated in a 
cylinder of radius 1 kpc centered on the bar. Vertical lines indicate 
the half-length of the gaseous component of the bar, determined 
by visual inspection. 



central region, even though only 1/4 of the available gas in 
the galaxy is located in the central region at that time. The 
remainder of the gas is mostly located at large radii, inside 
the spiral arms, where the gas density is significantly lower 
than in the central region. 

In the middle left and bottom left panels of Figure [5] 
we also plotted the oxygen mass Afo and iron mass Mpe 
inside the central region, after adjusting the scales to match 
the evolution of Mh at early time. Before t ~ 0.6 Gyr, the 
masses of hydrogen and oxygen increase roughly at the same 
specific rate M~^dM/dt, and as a result the [O/H] remain 
fairly constant (top right panel). The same is true for iron, 
up to t ^ 0.8 Gyr (see middle right panel). Afterward, Mq 
and Mpe increase faster than Mh, causing the increase in 
[O/H] and [Fe/H]. Mh reaches a peak at t = 0.85 Gyr while 
Mo and Mpe peak slightly later. After reaching the peak, 
all masses decrease, but AIh decreases faster than Mo and 
Afpo, and consequently [O/H] and [Fe/H] keep increasing. 

The plots of [O/H] and [Fe/H] show a similar behavior: 
there is an initial drop in [O/H] and [Fe/H] caused by the 
inflow of low-metallicity gas into the center. This is partly 
compensated by enrichment from Type H SNe, which pro- 
duce some iron and a lot more oxygen. At f = 1 Gyr, the 
SFR reaches its peak and starts decreasing. This reduces 
the production of oxygen by Type H SNe, and causes a 
change of slope in [O/H]. In the bottom right panel of Fig- 
ureO we plot [0/Fe], which gives the relative importance of 
Type II and Type la SNe. There is a clear change of regime 
at t ~ 1.4 Gyr, when enrichment by Type la SNe takes over. 
Note that chemical enrichment is not entirely caused by stars 
formed during the simulation. The stars present in the initial 
conditions contribute as well. Therefore, enrichment by both 
Type II and Type la Sne starts at t = 0, following the age 
distribution of the existing stars. Note also that Type la SNe 
have a delay of ab out 0.4 Gyr from the birth of their progen- 
itors (see Fig. 1 of lKobavashi et al...20 00). As a result, when 
the SFR increases, the rate of Type la SNe also increases 



© XXX RAS, MNRAS 000,[THT8] 



8 Hugo Martel et al. 



all different processes 




1 1 1 1 1 1 1 1 


inflow + outflow 


SF -(- winds-SNe 
1 , . , . 1 . , . , 1 . , . , 




time [Gyr] 

Figure 6. Rate of change of hydrogen mass (black hncs and left 
axes) and oxygen mass (red lines and right axes) inside the cen- 
tral region, versus time, for Run F. Top panel: solid lines show gas 
flowing into the central region (top lines) and out of the central re- 
gion (bottom lines). Dotted lines show gas added by stellar winds 
and SNe outflows (top lines) and gas removed by star formation 
(bottom lines) . Middle panel: sum of the quantities plotted in the 
top panel. Solid lines show the net flux across the boundary of 
the central region; dotted lines show the net effect of matter inter- 
change between the gas and stellar components. Bottom panel: 
sum of the quantities plotted in the middle panel, showing the 
net rate of change of mass, combining all processes. Dashed lines 
indicate a constant mass. 



but with a 0.4 Gyr delay relative to the rate of Type II SNe. 
After the SFR reaches its peak and starts decreasing, the 
iron production, mainly due to Type la SNe, can exceed the 
oxygen production with a delay (depending on how quickly 
the SFR decreases) . Enrichment by Type la SNe is more ex- 
tended in time, and as a result [Fe/H] keeps increasing and 
[0/Fe] keeps decreasing, at least over the timescale of the 
simulation. However, note that we see an initial increase of 
[0/Fe], because we set the initial [O/Fe] to 0. If the initial 
value of [O/Fe] was significantly higher, which would likely 
be the case in high-redshift galaxies, we would not expect to 
see an initial upturn of [O/Fe] due to the increase of SFR. 
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Figure 7. Top panel: [O/H] in the gas phase in the central region, 
versus time, for Run F. Bottom panel: rate of change of hydrogen 
mass (black lines and left axis) and oxygen mass (red lines and 
right axis) inside the central region, versus time. These panels 
are copied from Figures \5\ and [6] respectively. Vertical dashed 
lines identify various characteristic epochs. Numbers in top panel 
indicate values of [O/H]. 
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Figure 8. Top panel: [Fe/H] in the gas phase in the central region, 
versus time, for Run F. Bottom panel: rate of change of hydrogen 
mass (black lines and left axis) and iron mass (green lines and 
right axis) inside the central region, versus time. Vertical dashed 
lines identify various characteristic epochs. Numbers in top panel 
indicate values of [Fe/H]. 



3.2.2 Metallicity Evolution 

Four processes can modify the mass of a particular element 
in the gas phase inside the central region: (1) gas flowing into 
the central region, (2) gas flowing out of the central region, 
(3) gas being converted into stars, and (4) gas being released 
by stars in the form of stellar winds and SNe ejecta (we will 
refer to these two processes together as stellar outflows),. We 
now investigate the role of these various processes. The top 



panel of Figure |6] shows the rate of change of hydrogen mass 
Mh (black lines and left axes) and oxygen mass AIq (red 
lines and right axes) inside the central region. The middle 
panel of Figure[B]shows the combined effect of matter flowing 
across the boundary of the central region (solid lines) and 
matter exchange between the stellar and gas phase (dotted 
lines). The bottom panel of Figure [S] shows the combined 
effect of all processes. 

Between t — 0.2 Gyr and t = 1.1 Gyr, there are large 
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Figure 5. Time-evolution of quantities inside the central region, for Run F. Top left panel: Star formation rate (solid line); middle and 
bottom left panels; hydrogen mass (black lines and left axes), oxygen mass (red line and right axis in middle panel), and iron mass (green 
line and right axis in bottom panel); right panels: [O/H], [Fe/H], and [0/Fo] in the gas phase, as indicated (solid lines). The dotted lines 
on four of the panels indicate the corresponding global quantities (see Fig.fSj. 
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Figure 9. Star formation history for Run F. Each dot corre- 
sponds to a star formation event (a gas particle turning into a 
star particle). The figure shows the epoch of the event versus ra- 
dius. The red line shows the boundary of the central region. The 
blue line shows the half-length of the gaseous component of the 
bar, determined by visual inspection. The dots located on the 
right side of that curve correspond to stars formed in the disc, 
mostly inside the spiral arms. 



amounts of gas both moving m and moving out of the cen- 
tral region (solid lines in top panel of Fig. [6]) . This can be 
understood by examining the velocity field of the gas inside 
the bar (Fig. [2]). At f = 0.5 Gyr, the gas inside the bar moves 
along an elongated, elliptical orbit. This orbit intersects the 
edge of the central region, hence the gas is flowing across the 
boundary, both moving in and moving out. This elongated 
vortex contracts with time, and after t ^ 1.1 Gyr, it is en- 
tirely contained inside the central region. At this point, the 
flux of gas through the edge is small, and further enrichment 
is caused mostly by stars in the central region (bottom dot- 
ted lines in top panel of Fig. |6} . Because of the overall con- 
traction of the gas along the bar, there is a net inflow of gas 
into the central region (middle panel of Fig. [S}. The rate of 
net inflow is quite irregular because the gas tends to become 
clumpy, but overall most of the inflow takes place between 
0.2 Gyr and 1.1 Gyr. As the gas get progressively enriched 
by SNe and stellar winds, the flux of oxygen increases rela- 
tive to the flux of hydrogen, causing a net increase of [O/H] 
inside the central region. 

The net effect of star formation, SNe, and stellar winds 
is a decrease in the amount of gas, since some of the mat- 
ter remains trapped in stellar remnants or in stars whose 
lifetime exceeds the duration of the simulation. Because gas 
returned to the ISM by SNe and winds has been enriched, 
the net effect is again an increase in [O/H]. There is a delay 



© XXX RAS, MNRAS 000, [THm 



10 Hugo Martel et al. 



t - □ 25 Gyr 

t = 0,50 Gyr 

t = 0.75 Gyr 

t - 1.00 Gyr 




R [kpc] 

Figure 10. Star formation rate (top panel) and specific star for- 
mation rate (bottom panel) versus radius, for Run F. The differ- 
ent curves correspond to different epochs, as indicated. Vertical 
color lines indicate the half-length of the gaseous component of 
the bar. The dotted black line indicates the radius of the central 
region. 



between gas flowing into the central region and gas being 
removed by star formation, the latter becoming important 
only after ~ 0.7 Gyr. As a result, the total mass of hydrogen 
and oxygen in the central region increases with time until 
t ^ 1.0 Gyr, and then decreases afterward (bottom panel of 
Fig.©. 

In Figure [T] we combine the top right panel of Figure [5] 
with the middle panel of Figure (6] to show the effect of 
these various processes on the central metallicity. Between 
t = 0.3 Gyr and t = 0.7 Gyr, [O/H] increases from 0.12 
to 0.22, even though the effect of the central stars is neg- 
ligible. This early enrichment is entirely caused by metal- 
rich gas flowing into the central region. This gas was not 
enriched by central stars, and was therefore enriched in 
other regions, prior to falling into the central region. Af- 
ter t = 0.7 Gyr, central star formation becomes important, 
but until t = 0.96 Gyr, gas flow across the boundary remains 
the dominant process. During this period, [O/H] increases 
from 0.22 to 0.32. Star formation starts to dominate only af- 
ter t = 0.96 Gyr, and gas flow across the boundary becomes 
negligible after t — 1.18 Gyr. Overall, half of the enrichment 
(from [O/H] = 0.12 to [O/H] = 0.32) takes place before 
t ~ 1.0 Gyr, when net inflow of metal-rich gas was either 
the only process or the dominant process taking place. In 
Figure [8l we plot the corresponding results for iron. This re- 
sults are qualitatively similar to the ones shown in Figure [T] 
Gas inflow dominates until f ~ 1 Gyr. Afterward, the dom- 
inant process is net removal of gas by star formation and 
stellar outflows. 

These results imply that the central region does not 
evolve in isolation. About half of the metals that end up in 
the central region were formed somewhere else. To under- 



stand this, we need to examine the global history of star 
formation in the galaxy. We identified all star formation 
events in the simulation (that is, a gas particle turning into a 
star particle between two output steps). Figure |9] shows the 
epoch and the distance from the center where each event 
took place. The red line indicates the radius of the central 
region. The blue curve shows the extent of the gaseous com- 
ponent of the bar, and was calculated separately by visual 
inspection of the right-hand-side panels of Figure[T]and sim- 
ilar panels at different times. The blue curve forms a near- 
perfect envelope that encloses most star formation events. 
This indicates that star formation is not limited to the cen- 
tral region. Instead, stars are forming everywhere along the 
bar. At t — 0.5 Gyr, the gaseous component of the bar is 
about 6 Gyr in length, that is, three times the diameter 
of the central region. Most of the stars forming inside the 
bar at that time form outside of the central region. As the 
gaseous component of the bar contracts, a larger fraction of 
the stars form inside the central region, and at t ~ 1.1 Gyr, 
the gaseous component of the bar falls entirely inside the 
central region. This corresponds to the time when the next 
flux of metal-enriched gas across the boundary of the central 
region becomes negligible, as seen in Figures [7] and |8] 

We calculated the SFR profiles and specific SFR (SSFR) 
profiles at different times. The results are shown in Fig- 
ure [TO] The curves at 0.50 Gyr are rather fiat, indicating 
that star formation is taking place all along the bar, at 
roughly the same rate0 Afterward, the gaseous component 
of the bar contracts, and the SFR and SSFR become cen- 
trally concentrated. However, the fiuctuations in SSFR at 
late time remain fairly small. The SFR becomes centrally 
concentrated at late time mostly because the gas available 
for forming stars becomes itself centrally concentrated. 

Notice that star formation is not limited to the bar. All 
the dots on the right of the blue curve in Figure[9]correspond 
to stars forming in the disk, mostly inside the spiral arms. 
The apparent "void" between 1 and 5 kpc, which is also vis- 
ible in the top red curve in Figure [TOl is simply a geometric 
effect; if stars are forming uniformly over the surface of the 
disk, we expect the SFR to increase linearly with radius. 



3.3 Changing the Resolution 

Figure [TT] shows the evolution of the SFR, hydrogen mass, 
[O/H], [Fe/H], and [0/Fe] in the central region, for all runs. 
The blue and black lines show the results for Runs HR 
and F, respectively, which only differ in the number of par- 
ticles used. As we see, increasing the resolution does not 
lead to convergence. This is likely caused by the use of 
a fixed density threshold, nth, for star formation. As the 
resolution increases, the simulation is capable of resolving 
finer structures, and more gas particles are eligible to reach 
the threshold density more quickly. The finer structures 
and earlier starbursts could help developing more struc- 
tures in the gas, which could enhance the growth of the 
non-asymmetric structures of the stars, such as a bar and 
spiral arms. We should be able to adjust nth and the other 
parameters to match the lower-resolution run, but this is 



^ At t = 0.25 Gyr, the bar has not formed yet. 
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Figure 11. Time-evolution of quantities inside tlie central region. Top left panel: Star formation rate; middle left panel: hydrogen mass; 
right panels: [O/H], [Fe/H], and [O/Fe] in the gas phase, as indicated. Black lines: Run F; cyan lines: Run HR; green lines: Run NoBar; 
red lines: Run Feed. 



not our goal in this paper. This comparison run demon- 
strates that the timescale of the evolution of chemical and 
dynamical properties is sensitive to the simulation set-up, 
but the evolution paths are qualitatively similar. In future 
work, we will continuously improve our simulations by cal- 
ibrating the se parameters against th e various observational 
constraints l|Rahimi &: Kawatall2012l ). and explore also the 
initial conditions, such as the mass of the disk and the gas 
mass fraction. 

Still, the results for the Runs F and HR are qualita- 
tively similar. The SFR and the hydrogen mass A/h initially 
increase, reach a peak, and then decrease. The SFR peak 
is reached earlier in Run HR, at f = 0.7 Gyr instead of 
1.0 Gyr. Looking at the right panels in Figure 1111 we see 
interesting differences. In Run HR, [O/H] reaches a peak 
at t ~ 1.5 Gyr, and starts to decrease. The plot of [O/Fe] 
versus time shows that after t = 1 Gyr, Type la SNe and 
mass loss from the intermediate-mass stars start to domi- 
nate. The stellar mass ejecta do not produce much oxygen, 
but they produce a lot of hydrogen. Hence, the decrease in 
[O/H] is not caused by oxygen depletion, but rather by dilu- 
tion: [O/H] goes down because the amount of hydrogen goes 
up relative to the amount of oxygen. In Run F, mass loss 
from intermediate-mass stars takes over much later, and for 



that reason [O/H] keeps increasing until the end of the sim- 
ulation. Figure [121 shows the evolution of [O/H], Mn, and 
Mo in the central region, for Run HR. The vertical dashed 
lines indicate the same particular epochs as in Figure [T] star 
formation becoming significant, equal contribution from in- 
flow and stars, and inflow becoming insignificant. Compar- 
ing with Figure[7l we see that the period during which inflow 
+ outflow of gas dominates is much shorter, extending from 
t — 0.20 Gyr to t = 0.63 Gyr. Also, star formation and stel- 
lar outflows become important earlier, around t = 0.38 Gyr. 
But qualitatively, the results are the same for Runs F and 
HR: an early period dominated by gas flowing in and out of 
the central region, during which [O/H] increases from 0.16 
to 0.32, followed by a late period dominated by star forma- 
tion and stellar outflows, during which [O /H] increases to a 
maximum of 0.56 before falling. 



3.4 Barred versus Unbarred Galaxy 

The green lines in Figure [TT] show the results for Run No- 
Bar. In the absence of a bar, there is very little gas moving 
into the central region, and both the SFR and Mh increase 
very slowly, resulting is smaller changes in [O/H], [Fe/H], 
and [O/Fe] compared to the fiducial run. Interestingly, at 
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Figure 12. Same as Figure [7] for Run HR. 



Figure 14. Same as Figure Ffl for Run Feed. 




time [Gyr] 

Figure 13. Same as Figure [7] for Run NoBar. Notice the dif- 
ference in vortical scales between the bottom panels of the two 
figures. 



t = 2 Gyr, the mass of hydrogen in the central region is al- 
most the same in Runs F and NoBar, though it gets to that 
value in totally different ways. In Run F, the bar channels 
large quantities of gas into the central region, but most of the 
gas is converted into stars. In Run Nobar, the SFR is much 
lower. The SFR increases as Mh increases, but never gets 
high enough to result in a decrease in Mh. Metal enrichment 
is reduced compared to Runs F and HR, and interestingly 
[O/Fe] is nearly constant. Figure [T51 shows the evolution of 
[O/H], Mh, and Mo in the central region, for Run NoBar. 
The results are totally different than the ones for Run F 
(Fig. El . The rates are about an order of magnitude smaller, 
for all processes. We do not find an early period dominated 
by inflow -f outflow and a late period dominated by star for- 
mation -I- stellar outflows. Instead, all processes are taking 
place simultaneously, and all rates increase roughly linearly 
with time. Inflow + outflow dominate at all times (bottom 



panel), indicating that the increase in [O/H] and [Fe/H] seen 
in the central region is mostly caused by gas being enriched 
outside of the central region and moving in afterward. In 
Runs F and HR, the gaseous component of the bar con- 
tracts and eventually falls entirely inside the central region, 
shutting down inflow. This does not happens in the absence 
of a bar. Gas flows slowly and steadily into the central region 
all the way to the end of the simulation. 



3.5 Stronger Feedback 

The red lines in Figure [TT] show the results for Run Feed. 
The results are very similar to Run F. Increasing the stel- 
lar wind and SNe feedback by an order of magnitude has 
virtually no effect on the flow of the gas, indicating that 
this flow is still dominated by gravitational dynamics. The 
only real difference is the central SFR and central hydro- 
gen mass, which start to rise earlier, but eventually reach 
the same peak values at the same time. Having nearly the 
same SFR in Runs F and Feed likely results from competing 
effects. A larger feedback results in compression and reheat- 
ing of the ISM, the former process favoring star formation 
and the second process inhibiting it. Still, it is somehow sur- 
prising that such a large increase in feedback makes so little 
difference. Figure [Til shows the evolution of [O/H], Mh, and 
Mo in the central region, for Run Feed. Again, the results 
are similar to those for Run F. Gas inflow dominates up 
to t = 1 Gyr, and during that period [O/H] increases from 
0.10 to 0.35. Afterward, star formation and stellar outflows 
dominate, and [O/H] increases up to 0.47. 



3.6 The Effect of Inclination and Aperture 

So far, we have defined the central region as a cylinder of 
diameter 2 kpc centered on the rotation axis of the galaxy. 
Observationally, this corresponds to a disc galaxy seen face- 
on. In practice, face-on galaxies are fairly rare. In this sec- 
tion, we redefine the central region as a cylinder of diameter 
2 kpc, not centered on the rotation axis, but rather on the 
line-of-sight. We consider galaxies with a 45° inclination. 
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Figure 15. Time-evolution of quantities inside the central region, for Run F. Solid lines: 45° inclination; dotted lines: no inclination. 
Top left panel: Star formation rate; middle and bottom left panels: hydrogen mass (black lines and left axes), oxygen mass (red line and 
right axis in middle panel), and iron mass (green line and right a^cis in bottom panel); right panels: [O/H], [Fe/H], and [O/Fe] in the gas 
phase, as indicated. On these three panels, the solid and dotted curves differ by only a few percents, and are visually indistinguishable. 



This is large enough to potentially have an effect, while still 
small enough to allow for a visual identification of the bar. 

Figure [T^ shows the results for Run F. For comparison, 
we also plot as dotted lines the results with no inclination, 
copied from Figure [S] They are very similar. In particular, 
[O/H], [Fe/H], and [O/Fe] differ by less than 2%, making 
the curves indistinguishable. Most of the gas in the central 
region is concentrated in the inner 1 kpc, so rotating a region 
of size 2 kpc only produces a small edge effect. Therefore, our 
results are essentially insensitive to inclination. 

We have considered an aperture of 2 kpc, because it cor- 
responds to the typical size of the fibre in the SDSS. But of 
course the actual aperture is redshift-dependent. We now go 
back to Run F with no inclination (face-on galaxy), and con- 
sider the effect of varying the aperture. Figure [16] shows the 
evolution of [O/H], Mh, and Mo in the central region, for 
an aperture of diameter 1 kpc. Comparing with Figure [71 
we find that reducing the aperture does not qualitatively 
change the results. We still have an early phase when gas 
inflow dominates, and a later phase when star formation and 
stellar outflows dominate. However, everything is shifted to 
later times. The effect of star formation and stellar outflows 
become important &t t = 0.8 Gyr (compared to 0.7 Gyr for 
the 2 kpc aperture), the effects of gas inflow and stars are 
comparable at t = 1.4 Gyr (compared to 0.96 Gyr for the 
2 kpc aperture), and at the end of the simulation, gas inflow 



is still important. The gaseous component of the bar con- 
tracts with time, but it eventually reaches a minimum size 
which roughly corresponds to the aperture size of 1 kpc (see 
Fig. [9]). 

Figure [17] shows the results for an aperture of diameter 
5 kpc, which encloses most of the bar for the full duration 
of the simulation. This time everything is shifted to earlier 
times. The effects of star formation and stellar outflows be- 
come important a.t t = 0.4 Gyr, the effects of gas inffow and 
stars are comparable at t = 0.73 Gyr, and gas inflow becomes 
negligible sX t = 0.96 Gyr, when the gaseous component of 
the bar falls entirely inside the central region. Interestingly, 
the evolution of [O/H] is very insensitive to the choice of 
aperture (top panels in Figs. [161 and I17|l . even though the 
timing of the physical processes responsible for the evolution 
of [O/H] strongly depends on the choice of aperture. 



4 DISCUSSION 

The classic picture of gas flowing along the bar, produc- 
ing a central starburst, and enriching the gas in the central 
region, turns out to be too simplistic. Several physical pro- 
cesses complicate this picture. The gas falling toward the 
central region tends to form an elongated inner orbit inside 
the stellar bar. Until this orbit contracts sufficiently to be 
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Figure 16. Same as Figure for Run F and a central region of 
diameter 1 kpc. Dotted line in top panel shows the evolution of 
[O/H] for an aperture of 2 kpc. 
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Figure 17. Same as Figures [7] and 1161 for Run F and a cen- 
tral region of diameter 5 kpc. Dotted line in top panel shows the 
evolution of [O/H] for an aperture of 2 kpc. 



entirely contained inside the central region, this region does 
not evolve in isolation. The dominant process during this 
period is not star formation + stellar outflows in the central 
region, but rather exchange of gas between the central region 
and the rest of the galaxy. The amount of gas flowing in and 
out of the central region greatly exceeds the amount of gas 
involved in star formation and stellar outflows, even when 
the central SFR reaches its maximum value. This implies, vn 
particular, that metal-enrichment in the central region can- 
not he viewed as a mere consequence of star formation m 
the same region. Stars formed in the central region can end 
up enriching the gas outside of the central region, and vice- 
versa. 

For instance, consider a star particle forming in the cen- 
tral region. This star particle represents collectively an en- 



Figure 18. Histograms of the number of times gas and star par- 
ticles have crossed the boundary of the central region between 
t = Gyr and t = 2 Gyr, for Run F. In the leftmost bin (zero 
crossings), we excluded particles that were never inside the cen- 
tral region. The seesaw pattern results from the fact that many 
particles start outside the central region at t = Gyr and end up 
inside at t = 2 Gyr, which implies an odd number of crossings. 



semble of stars which have a range of initial masses, and 
therefore a range of lifetimes. The production of metals by 
these stars and their deposition into the ISM will therefore 
be extended in tim^f]. As the star particle moves along the 
length of the bar on an elongated orbit, it will deposit met- 
als over the entire length of the bar. The same is true for 
star particles forming in the bar but outside of the central 
region: they will deposit metals over the entire length of the 
bar, including in the central region. Furthermore, the gas it- 
self follows elliptical orbits, and this contributes to spreading 
metals over the length of the bar. 

To illustrate this, we calculated the number of times 
each particle crossed the boundary of the central region, 
in Run F. Figure [TS] shows histograms of the number of 
crossings, for star and gas particles. We excluded particles 
that were never inside the central region. Hence, the leftmost 
bins (zero-crossings) represent particles that were located in 
the central region in the initial conditions, and remained in 
the central region throughout the simulation. 

Star particles experience lots of crossings. The median 
value is 15 crossings. When a particle moves along an elon- 
gated orbit that intersects the central region, there are 4 
crossings per orbit. Hence, about half of the star particles 
complete between 1 and 4 full orbits. The histogram for gas 
particles is significantly diff'erent. While the length of stel- 
lar component of the bar remains roughly constant during 
the simulation, the gaseous component of the bar contracts 



^ Incidently, this explains why the abundance ratios vary much 
more smoothly than the SFR. 
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with time. As a result, most gas particles that cross the 
boundary of the central region cross it only once. Most of 
the gas located in the bar will essentially end up in the cen- 
tral region if we wait long enough. This implies that there 
is no correlation between central star formation and cen- 
tral metal-enrichment, but, in galaxies that are sufficiently 
evolved, there is a correlation between global star formation 
and central metal-enrichment, in the sense that stars form- 
ing all over the length of the bar produce and release metals 
which eventually end up in the central region. 

We focused on the central 2kpc because of its obser- 
vational relevance. It is interesting to see what an observer 
could infer about the evolutionary stage of the galaxy, sim- 
ply by observing the central region. All quantities plotted 
in Figure [5] (SFR, Mh, Mo, Mpe, [0/H], [Fe/H], [O/Fe]) 
vary non-monotonically with time. Hence, measuring one of 
these quantities could not determine the current evolution- 
ary stage of the galaxy. Only a combination of measurements 
could achieve that. To illustrate this, we plot the central 
value of [0/H] versus the central value of the SFR for Run 
F in Figure [T9l The line does not intersect itself, indicating 
that each point on the line corresponds to a specific epoch 
in the evolution of the galaxy. We have identified seven par- 
ticular epochs, and plotted the stellar and gas distributions 
for each one. 

• t = 0.0 Gyr. Beginning of the simulation. The stellar 
and gas distributions are axisymmetric, the SFR is very low, 
and [0/H] has the initial value. 

• t = 0.5 Gyr. This is the epoch when the bar is the 
strongest and there is plenty of gas. The central SFR is still 
very low, and the dominant process responsible for metal- 
enrichment of the central region is gas exchange between the 
central region and the rest of the bar (Fig. |6] top and middle 
panels; Fig. [T] bottom panel). 

• t = 0.7 Gyr. The flux of gas in and out of the central 
region gets close to its maximum value, and star formation 
-I- stellar outflows start to become important (Fig. [71 bottom 
panel) . 

• t — 0.96 Gyr. The SFR and gas mass in the central re- 
gion reach their maximum values (Fig. [5l left panels). The 
contribution of gas inflow and star formation -I- stellar out- 
flows to metal-enrichment of the central region are compa- 
rable (Fig.[71 bottom panel). The orbit of the gas in the bar 
will soon be enclosed by the 2 kpc aperture (Fig. (2] bottom 
left panel). 

• t = 1.2 Gyr. The orbit of the gas has fallen inside the 
central region, and the flux of gas across the boundary of the 
central region becomes insigniflcant. (Fig. [T] bottom panel). 
From that point on, the central region evolves like a closed- 
box system, and star formation -I- stellar outflows are solely 
responsible for further metal-enrichment of the central re- 
gion. 

• t = 1.5 Gyr. [O/Fe] reaches its peak and starts decreas- 
ing (Fig. [5] bottom right panel), indicating a transition from 
enrichment dominated by Type II SNe to enrichment dom- 
inated by Type la SNe. 

• t — 2.0 Gyr. End of the simulation. The SFR steadily 
decreases as the available gas in the central region is depleted 
(Fig. [5] left panels). [O/H] and [Fe/H] both increase at a 
constant rate (Fig. O top right and middle right panels). 



5 SUMMARY AND CONCLUSION 

Our objective was to investigate the connection between the 
central SFR and central metallicity in barred disc galaxies. 
In this initial study, we considered the particular case of a 
disc galaxy with a baryonic mass of 5 x 10^°Mq and initial 
stellar fraction of 80%, in which a bar forms by bar instabil- 
ity. We performed a series of four simulations. We analyzed 
these simulations, focusing on the evolution of the SFR and 
metallicity in the central region of the galaxy, defined as the 
central kiloparsec (or 2 kpc in diameter), the approximate 
size of the SDSS fibre. 

In our first simulation, a strong bar of length 9 kpc com- 
posed of stars and gas forms around t ~ 0.5 Gyr. The stars 
and gas inside the bar do not lose all their angular momen- 
tum. This prevents them from falling directly into the cen- 
ter. Instead, they form elongated orbits that go through the 
center of the galaxy. The orbits tend to contract with time, 
but the effect is much more important for the gas than the 
stars. By the end of the simulation, at t — 2 Gyr, the stellar 
components of the bar still has a length of order 7 kpc, while 
the gaseous component of the bar has contracted down to a 
length of less than 1 kpc. 

Before t = 0.5 Gyr, when the bar has not yet formed, 
star formation is taking place mostly inside a radius of 5 kpc, 
where the gas is the densest. At t = 0.5 Gyr, when the bar 
forms, 69% of the star formation takes place inside the bar, 
the remainder taking place inside the spiral arms. At this 
point, the gaseous component of the bar has a length of 
about 7 kpc, much larger than the size of the central region. 
Hence, most stars forming m the bar do not form in the cen- 
tral region. But as the gaseous component of the bar con- 
tracts, a larger fraction of that gas is located in the central 
region, and the SFR in that region increases. By t = 1.1 Gyr, 
that gas is enclosed entirely inside the central region. There 
is no more gas fiowing in, and the central SFR decreases as 
less gas is available to form new stars. 

We tracked the evolution of the metal abundances in 
the gas phase, in the central region. We start the simulation 
with an initial metallicity gradient. This gradient is main- 
tained in the early stages of the simulation, when gas moves 
along circular orbits. Once the bar forms, gas starts moving 
radially. The central metallicity initially decreases, and then 
steadily increase until the end of the simulation. Naively one 
could think that this increase in central metallicity is caused 
by outflows from central stars. We decided to investigate this 
issue in more detail, by calculating separately the contribu- 
tion of the various process that can change the hydrogen 
mass, oxygen mass, and iron mass located in the central 
region (Figs. [S] [T] andlH}. The surprising result is that the 
central metallicity increases significantly before the contribu- 
tion of central stars becomes significant. Star formation and 
metal enrichment is taking place along the entire length of 
the bar. Because initially the central region is much smaller 
than the bar, most of the metals are produced outside the 
central region. As the gas moves along the bar on ellipti- 
cal orbits, it moves in and out of the central region. But as 
these elliptical orbits contract with time, there is a net flux 
of gas moving into the central region. This gas was predom- 
inantly enriched by stars located outside the central region, 
and early-on, this net influx of metal-enriched gas into the 
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Figure 19. [O/H] in the gas phase versus SFR, in the central region, for Run F. Insets show configuration of star particles (blacls 
dots), gas particles (cyan dots), and the boundary of the central region (orange circles) at various times, as indicated. Each inset is 
24 kpc X 24 Icpc. For clarity, only 1/6 of the particles are plotted. 



central region dominates over the metals produced by the 
central stars. 

Overall, about half of the increase in [O/H] and [Fe/H] 
in central region is caused by in-situ stellar enrichment. The 
other half is caused by stars located outside the central re- 
gion, which produce metals that are carried into the cen- 
tral region by large-scale gas flows along the bar. The main 
conclusion is that there is no direct connection between cen- 
tral SFR and central metallicity, but there is a connection 
between global SFR and late-time, central metallicity: stars 



forming over the entire length of the bar produce metals that 
will eventually end up in the central region as the gaseous 
component of the gas contracts. 

Some recent IFU observations of barred galaxies appear 
to support the conclusions of our simulations that star for- 
mation happens along the leng th of the bar l| Cant in et al.l 
I2OI0I : [Robert et al.ll201ll . |2013| '1. Histograms of the mass of 
the stellar populations versus age for two barred galaxies, 
NGC4900 and NGC5430, reveal that the stellar metallicity 
in the central region does not simply increase gradually with 
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time. Although the stellar metallicity may be difFerent than 
the gas metallicity, it suggests that the gas used to form the 
stars in the central region may get enriched and diluted by 
an outside contribution (like gas flowing along the bar). Or 
the mixing process for the gas (between the time massive 
stars release enriched elements and low mass stars release 
less enriched elements) involves a complicated timing with 
the formation of new generation of stars. 

We performed three more simulations. In one case, we 
increased the resolution by a factor of 3, and in another case, 
we increased the amount of stellar wind and SNe feedback 
by a factor of 10. In both cases, we obtained results that are 
quantitatively different, but qualitatively similar to our ini- 
tial run. In particular, we still identify an early phase when 
net gas inflow into the central region is primarily responsi- 
ble for increasing the central metallicity, and a later phase 
when metal-enrichment by central stars dominate. We also 
performed a simulation in which no bar forms. In that sim- 
ulation, there is very little gas flowing across the boundary 
of the central region, and very few stars forming in the cen- 
tral region. But interestingly, gas inflow is still the dominant 
process responsible of enriching the gas in the central region. 

In this paper, we focused on a particular set of initial 
conditions, in order to investigate the relationships between 
the various physical processes driving the evolution of the 
galaxy. In future work, we will vary the initial conditions, 
specifically the mass of the galaxy and the strength of the 
bar, with the goal of explaining the origin of the various 
trends revealed by studies of barred galaxies in the SDSS. 
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